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4 Abstract. Numerical studies of the interplanetary "multiple magnetic 

5 clouds (Multi-MC)" are performed by a 2.5-dimensional ideal magnetohy- 

6 drodynamic (MHD) model in the heliospheric meridional plane. Both slow 

7 MCI and fast MC2 are initially emerged along the heliospheric equator, one 

8 after another with different time interval. The coupling of two MCs could 

9 be considered as the comprehensive interaction between two systems, each 

10 comprising of an MC body and its driven shock. The MC2-driven shock and 

11 MC2 body are successively involved into interaction with MCI body. The 

12 momentum is transferred from MC2 to MCI. After the passage of MC2-driven 

13 shock front, magnetic field lines in MCI medium previously compressed by 

14 MC2-driven shock are prevented from being restored by the MC2 body push- 
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15 ing. MCI body undergoes the most violent compression from the ambient 

16 solar wind ahead, continuous penetration of MC2-driven shock through MCI 

17 body, and persistent pushing of MC2 body at MCI tail boundary. As the 

18 evolution proceeds, the MCI body suffers from larger and larger compres- 

19 sion, and its original vulnerable magnetic elasticity becomes stiffer and stiffer. 

20 So there exists a maximum compressibility of Multi-MC when the accumu- 

21 lated elasticity can balance the external compression. This cutoff limit of com- 

22 pressibility mainly decides the maximally available geoeffectiveness of Multi- 

23 MC, because the geoeffectiveness enhancement of MCs interacting is ascribed 

24 to the compression. Particularly, the greatest geoeffectiveness is excited among 

25 all combinations of each MC helicity, if magnetic field lines in the interact- 

26 ing region of Multi-MC are all southward. Multi-MC completes its final evo- 

27 lutionary stage when the MC2-driven shock is merged with MCl-driven shock 

28 into a stronger compound shock. With respect to Multi-MC geoeffectiveness, 

29 the evolution stage is a dominant factor, whereas the collision intensity is 

30 a subordinate one. The magnetic elasticity, magnetic helicity of each MC, 

31 and compression between each other are the key physical factors for the for- 

32 mation, propagation, evolution, and resulting geoeffectiveness of interplan- 

33 etary Multi-MC. 
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1. Introduction 

34 Space weather refers to the conditions on the Sun and in the solar wind, magnetosphere, 

35 ionosphere, and thermosphere that can influence the performance and reliability of space- 

36 borne and ground-based technological systems or can endanger human life or health, as 

37 defined in US National Space Weather Program Implementation Plan. A seamless fore- 

38 casting system for Space weather lies on the comprehensive and in-depth understanding 

39 of the Sun-Earth system. The never-stopping tremendous efforts have been made by hu- 

40 mankind since the space age of the 1950s. A great deal of the sophisticated observations 

41 beyond the Earth are now provided, with the launching of various spacecraft into deep 

42 space, such as Yohkoh, Geotail, Wind, SOHO, Ulysses, ACE, TRACE in the 1990s, and 

43 Cluster, RHESSI, SMEI, DS, Hinode (Solar B), STEREO in the 21st century. These 

44 spacecraft missions construct an indispensable backbone in the establishment of space 

45 weather prediction system. Meanwhile, many models have been or are being developed 

46 and applied to space weather forecasting by utilizing most measurements of the above 

47 spacecraft, such as (1) HAF (Hakamada-Akasofu-Fry) [Fry et al, 2001, 2005]; (2) STOA 

48 (Shock Time of Arrival) [Smart and Shea, 1985]; (3) ISPM (Interplanetary Shock Propaga- 

49 tion Model) [Smith and Dryer, 1990]; (4) an ensemble of HAF, STOA and ISPM models 

50 [Dryer et al, 2001, 2004; McKenna-Lawlor et al, 2006]; (5) SPM (Shock Propagation 

51 Model) [Feng and Zhao, 2006]; (6) SWMF (Space Weather Modeling Framework) [Toth 

52 et al, 2005]; (7) HHMS (Hybrid Heliospheric Modeling System) [Detman et al, 2006]; 

53 (8) a data-driven Magnetohydrodynamic (MHD) model of the University of Alabama in 

54 Huntsville [Wu et al, 2005a, 2006a]; (9) a 3D regional combination MHD model with in- 
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55 puts of the source surface self-consistent structure based on the observations of the solar 

56 magnetic field and K-coronal brightness [Shen et ai, 2007]; (10) A merging model of SAIC 

57 MAS and ENLIL Heliospheric MHD Model [Odstrcil et ai, 2004b]; (11) an HAF + 3-D 

58 MHD model [Wu et ai, 2005c, 2006c, 2007b, c], and so on. However, great challenges are 

59 still faced to improve the prediction performance of space weather, as human civilization 

60 is relying more and more on space environment [Baker, 2002; Fisher, 2004]. 

61 The interplanetary (IP) space is a pivot node of the solar-terrestrial transport chain. 

62 Solar transients, e.g., shocks and coronal mass ejections (CMEs), propagate in it, interact 

63 with it, and cause many consequences in the geo-space. Magnetic clouds (MCs) are 

64 an important subset of interplanetary CMEs (ICMEs), occupying the fraction of nearly 

65 ~ 100% (though with low statistics) at solar minimum and ~ 15% at solar maximum 

66 [Richardson and Cane, 2004, 2005], and have significant geoeffectiveness [Tsurutani et ai, 

67 1988; Gosling et ai, 1991; Gonzalez et ai, 1999; Wu and Lepping, 2002a, b; Wu et ai, 

68 2003, 2006b; Huttunen et ai, 2005]. The current intense study of MCs could be traced 

69 back to the pioneer work by Burlaga et al. [1981], who firstly defined an MC with three 

70 distinct characteristics of enhanced magnetic field strength, smooth rotation of magnetic 

71 field vector, and low proton temperature, and described it as a flux rope structure. An 

72 MC is widely thought to be the IP manifestation of a magnetic flux rope in the solar 

73 corona, which loses equilibrium and then escapes from the solar atmosphere into the IP 

74 space [Forbes et al., 2006], with its both ends still connecting to the solar surface [Larson 

75 et al, 1997]. 

76 It is very likely for solar transients to interact with each other on their way to the Earth, 

77 especially at solar maximum when the daily occurrence rate of CMEs is about 4.3 in aver- 
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78 age on basis of the SOHO/Lasco CME catalogue (http://cdaw.gsfc.nasa.gov/CME_list). 

79 Some IP complicated structures were reported, such as complex ejecta [Burlaga et al, 

80 2002], multiple MCs (Multi-MC) [Wang et al, 2002, 2003a], shock-penetrated MCs [Wang 

81 et al., 2003b; Berdichevsky et ai, 2005; Collier et al., 2007], non-pressure-balanced "MC 

82 boundary layer" associated with magnetic reconnection [Wei et al., 2003a, b, 2006], ICMEs 

83 compressed by a following high-speed stream [Dal Lago et al., 2006], multiple shock in- 

84 teractions [Wu et al, 2005d, 2006d, 2007a]. However, all space-borne instruments, except 

85 the heliospheric imagers onboard SMEI and STEREO, observe either the solar atmo- 

86 sphere within 30 solar radii by remote sensing, or the in-situ space by local detecting, or 

87 both. Thus, numerical simulations are necessary to understand the whole IP dynamics. 

88 Below is an incomplete list of numerical studies of dynamical processes of CMEs/MCs 

89 and complex structures in the IP medium mentioned before: an individual CME/MC 

90 [Vandas et al, 1995, 1996, 2002; Groth et al, 2000; Schmidt and Cargill, 2003; Odstrcil 

91 et al, 2003, 2004a, 2005; Manchester et al, 2004; Wu et al, 2005b], the interaction of a 

92 shock wave with an MC [Vandas et al, 1997; Xiong et al, 2006a, b], the interaction of 

93 multiple shocks [Wu et al, 2004a, b, 2005d, 2006d, 2007a], and the interaction of multiple 

94 ejecta [Gonzalez- Esparza et al, 2004; Gonzalez- Esparza, 2005; Lugaz et al, 2005; Xiong 

95 et al, 2005; Wang et al, 2004, 2005b; Wu et al, 2006c, 2007c; Hayashi et al, 2006]. 

96 Therein, Wu et al [2005d, 2006d, 2007a] performed a 1.5-D MHD model to simulate the 

97 famous Halloween 2003 epoch, in which eruption time of solar flares was used as input 

98 timing for solar disturbances to study the shock-shock interaction (and overtaking) and 

99 the matching of shock arrival time at 1 AU with observations (ACE). In addition, Wu 

100 et al. [2006c, 2007c] performed 3-D global simulations by combining two simulation models 
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101 (HAF + 3-D MHD) to study the interacting and overtaking of two ICMEs. These obser- 

102 vation and simulation efforts do advance our understanding of solar-terrestrial physics. 

103 The Multi-MCs have already been verified by observations to be an important IP origin 

104 for the great geomagnetic storms [Wang et al, 2002, 2003a; Xue et ai, 2005; Farrugia 

105 et al., 2006; Xie et al., 2006; Zhang et al., 2007]. Particularly, for the 8 extremely large 

106 geomagnetic storms with Dst < —200 nT during the year 2000 ~ 2001, 2 of them were 

107 caused by Multi-MCs and one caused by shock-MC interacting structure [Xue et al., 2005]. 

108 Most recently, via summarizing the efforts of the NASA Living With a Star (LWS) Co- 

109 ordinated Data Analysis Workshop (CDAW) held at George Mason University, in March 
no 2005, Zhang et al. [2007] proposed that 24 out of 88 (27%) major geomagnetic storms 

111 with Dst < —100 nT from the year 1996 to 2005 were produced by multiple interacting 

112 ICMEs arising from multiple halo CMEs launched from the Sun in a short period. So the 

113 Multi-MC plays a notable role in producing large geomagnetic storms. There are two pos- 

114 sible conditions for double-MC formation [Wang et al., 2004]: (1) The speed of following 

115 MC should be faster than that of preceding MC; (2) The separation between the eruption 

116 of two MCs should be moderate (about 12 hours based on statistics of observed events). 

117 Evolutionary signatures of ICMEs interacting are found from spacecraft observations, i.e., 

118 heating of the plasma, acceleration/deceleration of the leading/trailing ejecta, compressed 

119 field and plasma in the leading ejecta, possible disappearance of shocks, and strengthening 

120 of the shock driven by the accelerated ejecta [Farrugia and Berdichevsky, 2004]. Previous 

121 simulations of interaction between two magnetic flux ropes in the IP space [Lugaz et al., 

122 2005; Wang et al., 2005b], the solar corona [Schmidt and Cargill, 2004; Wang et al., 2005a; 

123 Lugaz et al., 2007], and a local homogeneous medium background [Odstrcil et al., 2003] 
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124 only address a few typical cases in the dynamical aspect. Here a comprehensive study 

125 of many cases of MCs interacting under various conditions is carried out for better un- 

126 derstanding of both dynamics and ensuing geoeffectiveness. The interaction between two 

127 systems, each comprising of an MC and its driven shock, could be considered in some 

128 senses as a generalization of our recent studies of MC-shock interaction [Xiong et al, 

129 2006a, b]. Thus we address the following two issues naturally: (1) What is the role of the 

130 following MC body in Multi-MC evolution in comparison with our previous studies [Xiong 

131 et al, 2006a, b] of MC-shock interaction? (2) At what evolutionary stage a Multi-MC at 

132 1 AU reaches the maximum geoeffectiveness? The above answers are explored by a 2.5-D 

133 numerical model in ideal MHD process. 

134 The force-free magnetic flux rope models have been proven to be very valuable to 

135 interpret in-situ observations of MCs [e.g., Lundquist, 1950; Burlaga, 1988; Farrugia 

136 et al, 1993; Chen, 1996; Owens et al, 2006]. Particularly, Lundquist model [Lundquist, 

137 1950] is adopted in our model to describe the magnetic field configuration of an MC, as 

138 widely applied in the space science literature [e.g., Vandas et al., 1995, 1996; Wang et al., 

139 2002, 2003d, 2005b; Xiong et al, 2006a, b]. A following fast MC overtaking and inter- 

140 acting a preceding slow one in the IP space could result in a Multi-MC structure [ Wang 

141 et al., 2002, 2003a]. In order to explore the basic physics process of Multi-MC, we make 

142 the following assumptions to simplify the complex circumstance of double-MC structure 

143 in the numerical MHD simulation: (1) two MCs' axes parallel or anti-parallel with each 

144 other; (2) their axes are both within the ecliptic plane and perpendicular to the Sun-Earth 

145 line; (3) each MC is symmetric in the azimuth direction of the heliosphere, and consid- 

146 ered as an ideal loop encompassing the Sun; (4) magnetic reconnection does not exist in 



DRAFT 



April 5, 2009, 4:06aiii 



DRAFT 



XIONG ET AL.: MC-MC INTERACTION AND ITS GEOEFFECTIVENESS X - 9 

147 double MC interacting; (5) both MCs have the same size, mass, magnetic field strength, 

148 and plasma (3. Thus, two MCs in our model only differ in magnetic helicity sign and 

149 initial radial lift-off speed Vmc- A parametric study of Hmc and Vmc is focused in our model 

150 for the very specialized Multi-MC structure. Since the two MCs are very alike except Hmc 

151 and Vmc, they could be, to some extent, considered to be identical. MCI and MC2 are 

152 respectively used to label the two MCs launched from the Sun, one after another. Because 

153 an MC boundary is a self-enclosed magnetic surface, and two MCs' magnetic field lines 

154 would not blend under the condition of the strictly ideal MHD process, the sub-structures 

155 of double MCs, corresponding to the previously separated MCI and MC2 before collision, 

156 could be easily differentiated, and accordingly named as sub-MCl and sub-MC2. 

157 The goal of the present work is to conduct a systematic investigation of Multi-MC in 

158 the IP space. We give a brief description of the numerical MHD model in Section 2, de- 

159 scribe the dynamical behavior of MC-MC interaction in Section 3, discuss the consequent 

160 geoeffectiveness in Section 4, analyze the compressibility of MC-MC collision in Section 

161 5, and summarize the paper in Section 6. 

2. Numerical MHD Model 

162 The Multi-MC simulation is accommodated by a few slight modifications from our 

163 previous numerical model for MC-shock interaction [Xiong et ai, 2006a, b]. These mod- 

164 ifications are as follows, (1) The top boundary of simulated domain is extended from 

165 300 to 400 Rs] (2) The following shock is replaced by a following MC; (3) The initial 

166 speed Vmc, emergence time tmc, and magnetic helicity Hmc out of all input parameters 

167 for each sub-MC of Multi-MC are independently selected to make various combinations 

168 for parametric study shown in Table 2. First, the propagation through the IP space is 
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169 modeled by numerical simulation. Then, the geomagnetic storm excited by the solar wind- 

170 magnetosphere-ionosphere coupling is approximated by an empirical formula of Burton 
^^df^^ ~ Qi^) ~ ^^^^ [Burton et al, 1975]. Here the coupling function Q = Vr- Mm{Bz, 0) 

172 and the diffusion time scale r = 8 hours, with the radial solar wind speed and south-north 

173 magnetic field component respectively denoted by and B^. Burton model [Burton et al, 

174 1975] for geomagnetic disturbance has been analyzed and validated [Wang et al, 2003d; 

175 Wang, 2003e], and applied in Dst evaluation [Wang et al, 2003d; Wang, 2003e; Xiong 

176 et al., 2006a, b]. Thus the physical process of cause-effect transport chain for solar dis- 

177 turbances is fully described in our model. Moreover, the MC2-driven shock in all of our 

178 simulation cases is faster than the local magnetosonic speed all the way, and strong enough 

179 so that it would not be dissipated in the low j3 MCI medium [Xiong et al, 2006a, b]. 

3. MC1-MC2 Interaction 

ISO All 48 simulation cases of MC1-MC2 coupling are assorted into 4 groups in Table 2, 

181 with 18 cases of individual MC in 2 groups of Table 1 for comparison. Here, IM, EID, 

182 CID respectively stand for "Indiviual MC", "Eruption Interval Dependence", "Collision 

183 Intensity Dependence" , with the subscripts 1 and 2 denoting the sign of magnetic helicity. 

184 Case Ci is shared by Groups EIDi and CIDi, and Case C2 by Groups EID2 and CID2. 

185 In our simulation, an MC with southward/northward magnetic field in its rear half is 

186 defined to have positive/negative helicity. Both MCs are associated with positive helicities 

187 in Groups EIDi and CIDi, meanwhile MCI and MC2 are respectively associated with 

188 positive and negative helicities in Groups EID2 and CID2. 

189 The numerical simulation is performed in the ideal MHD process. The artificial nu- 

190 merical magnetic reconnection between MCs is strictly ruled out by a specific numerical 
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191 technique [c.f. Hu et ai, 2003]. Thus the dynamics in Groups IMi, EIDi, and CIDi is 

192 nearly the same as that in Groups IM2, EID2, and CID2 respectively, whereas the geoef- 

193 fectiveness is highly different due to the reversed north and south magnetic components 

194 within the cloud with opposite helicity. Moreover, by changing Dt [Dt = tmc2 — tmci, 
"5 tmci = hour), the initiation delay between a preceding MC of 400 km/s and a following 
i9e MC of 600 km/s in Groups EIDi and EID2, the Multi-MC formed by the MCI and MC2 

197 may reach different evolutionary stages on its arrival at 1 AU. Therefore the eruption 

198 interval dependence for MC1-MC2 interaction is easily discriminated by a comparative 

199 study. Similarly, collision intensity dependence is also explored by a parametric study of 
2™ Vmc2 from 450 to 1200 km/s in Groups CIDi and CID2. Meanwhile the full interaction 

201 between sub-clouds within 1 AU to maximally highlight collision effect is guaranteed by 

202 tmc2 = 12.2 hours in Groups CIDi and CID2. Furthermore, an individual MC with its 

203 speed from 400 to 1200 km/s in Groups IMi and IM2 supplements indispensably to other 

204 Groups for the study of coupling effect of two MCs. Cases Bi and B2 with tmc2 = 30.1 

205 hours, Ci and C2 with tmc2 = 12.2 hours, are typical examples of Multi-MC in the early 

206 and late evolutionary stages respectively, which are addressed below in details. 

3.1. Case Bi 

207 In Case Bi, we discuss the results of MC1-MC2 interaction for eruption speed Vmci = 

208 400 km/s, Vmc2 = 600 km/s, and initiation delay tmc2 = 30.1 hours. Figure 1 shows 

209 the successive behavior of MC1-MC2 interaction of Case Bi. The magnetic field lines, 

210 among which two are enclosed white solid lines marking the boundaries of MCI and MC2 

211 respectively, are superimposed on each color-filled contour image, and two radial profiles, 

212 one through the equator (noted by Lat. = 0°), the other through 4.5° southward (white 
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213 dashed lines in the images, noted by Lat. = 4.5°S), are plotted below. One can read 

214 the global vision from the images and local details from the profiles simultaneously for 

215 the propagation and evolution of Multi-MC. For better highlighting the local disturbance, 

216 Figures l(a)-(c) show the magnitude B of magnetic field from which the initial value 

217 B\t=o is deducted. Two identical MCs are successively injected into the IP space with 

218 different initial eruption speed. As long as the fast following MC2 lags behind the slow 

219 preceding MCI, each of them behaves as an individual event, and satisfies the criteria of a 

220 single MC. Because the MC-driven shock and incidental shock [Xiong et al, 2006a] both 

221 propagate along the heliospheric current sheet (HCS) in the IP medium, their inherent 

222 traits are identically characterized by a concave-outward morphology with the position 

223 of the strongest intensity being roughly 4.5° away from the HCS. MC2-driven shock just 

224 approaches MCI body tail at 46.5 hours, as seen from Figure 1(d). Across this shock 

225 front, radial speed Vr increases abruptly from 440 km/s at MCI tail to 670 km/s at MC2 

226 head. From then on, MC2 and MCI will directly collide to form a special IP complex 

227 named Multi-MC by Wang et al. [2002, 2003a], and their evolution will be coupled with 

228 each other. Consequently, the characteristic parameters of each sub-MC would change 

229 drastically due to the non-linear interaction. At 56.1 hours, MC2-driven shock front 

230 has already entered MCI body across which radial speed Vr abruptly jumps from 445 

231 to 620 km/s, but MC2 body is still unable to catch up with MCI tail (Figure 1(e)) 

232 because of tmc2 = 30.1 hours. The dynamic response of Multi-MC at this snapshot is 

233 merely ascribed to the interaction between MC2-driven shock and MCI body. So the 

234 preceding MCI behavior in Figures 1(b), (e), and (h) are similar to its counterpart of 

235 MC-shock interaction in essence [c.f. Figures 3(c), (f), and (i) in Xiong et al., 2006a]. 
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236 Large compression within MCI medium downstream of MC2-driven shock front is very 

237 pronounced from an abnormal local spike-like structure of Cf along Lat. = 4.5°S, as shown 

238 in Figure 1(h). The orientation of magnetic field lines is also rotated in MCI medium swept 

239 by the shock front. As the shock continuously advances into MCI body, the morphology 

240 of MCI rear part is transformed from an original rough semi-circle (Figure 1(d)) to a V- 

241 shape with a wide open mouth (Figure 1(f)). Moreover MC2 body has already contacted 

242 MCI tail at the bottom of so-called V shape along the equator at 80.7 hours, when the 

243 MC2-driven shock cannibalizes the rear half of MCI body (Figure 1(f)). Since then, MC2 

244 body is directly involved into interaction with MCI body. The Multi-MC evolution has 

245 reached a new critical stage, for MCI will undergoes the most violent compression from 

246 the ambient solar wind ahead, continuous penetration of MC2-driven shock through MCI 

247 body, and persistent pushing of MC2 body at MCI tail boundary. In Figure 1(f), nearly 

248 constant speed in MCI rear half and large speed difference with 80 km/s across MCI 

249 rear boundary along the equator imply continuous strike of high-speed MC2 body upon 

250 preceding MCI body. Besides, the interplanetary magnetic field (IMF) within Multi-MC 

251 envelope is highly bending just behind MC2-driven shock front (Figures 1(c), (f), and (i)), 

252 as a result of rotation across the shock front and draping around either sub-cloud surface. 

253 The in-situ observation along Lat. = 4.5°S by a hypothetical spacecraft at Lagrangian 

254 point (LI) is illustrated in Figure 2. With each sub-MC boundary identified as dashed 

255 lines, the MCI duration of 18 hours is much less than MC2 duration of 26 hours due to the 

256 compression in MCI rear half accompanying with MC2-driven shock advancing. The MC2 

257 "senses" the existence of preceding MCI, though its response is much less sensitive. The 

258 location of maximum bulk flow speed Vr in MC2 body is shifted by 6 hours later (Figure 
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259 2(c)) from MC2 head [c.f. Figure 2 in Xiong et al, 2006a], between which magnitude B 

260 is obviously enhanced (Figure 2(a)). The dawn-dusk electric field VB^ is calculated by 

261 the product of tv and Bg in the spherical geometry of this simulation. Beginning from 

262 74 hours, VBz, negative in MCI rear half, positive in MC2 front half, and negative again 

263 in MC2 rear half (Figure 2(d)), is respectively responsible for Dst dropping from nT 

264 at 74 hours to —140 nT at 82 hours, recovering from —140 nT at 82 hours to —25 nT at 

265 97 hours, and dropping again from —25 nT at 97 hours to —75 nT at 114 hours (Figure 

266 2(e)). Owing to compression of southward magnetic component Bg (Bg = Mm{Bg,0)) 

267 within MCI rear part, the first Dst dip with —140 nT is much lower than the second one 

268 with —75 nT for geomagnetic storm. Particularly, the two Dst dips are separated by only 

269 32 hours, because the geoeffectiveness of two IP triggers (MCI and MC2) are superposed 

270 together. The idea of a two-ejecta event associated with a two-step geomagnetic storm 

271 was recently proposed and verified by Farrugia et al. [2006] on basis of observation. Hence 

272 the association of two Dst dips lies in the MC1-MC2 interaction. 

3.2. Case Ci 

273 In order to realize the fully interaction between MCI and MC2 before their arrival at 

274 LI, tmc2; the emergence time of MC2, is scheduled earlier to be 12.2 hours with both 

275 MCs having the same speeds of Case Bi. Only the evolution of Vr is given in Figure 3 

276 to visualize multi-cloud structure. Comparing to that in Figures 1(c), (f), and (i), the 

277 so-called "V-shape" morphology of MCI rear half becomes very fiat under the pounding 

278 of very high-speed MC2 body at 19.5 hours as Multi-MC evolution proceeds, as shown in 

279 Figure 3(a). As a result, contact position between MCI and MC2 body is extended from 

280 one single point at the HCS (Figure 1(f)) to a straight line between Lat. = 4.5°S and 
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281 4.5°N (Figure 3(a)). The MCl's magnetic elasticity seems to be too vulnerable to resist the 

282 violent collision from MC2 body. The collision efficiently transfers the radial momentum 

283 from the fast following MC2 to the slow preceding MCI. It results in monotonically 

284 decreasing Vr from the head to tail of Multi-MC at 53.3 hours, resembling a single MC, 

285 as seen in Figure 3(c). Besides, MC2 morphology turns from a radial-extent-elongated 

286 ellipse (Figure 3(a)) to an angular-extent-elongated one (Figure 3(c)) due to the blocking 

287 of MCI body ahead. MC2 body is also compressed radially to some extent. Certainly, 

288 the compression of MC2 body is much less than that of MCI body. Moreover, MC2- 

289 driven shock ultimately penetrates the MCI body (Figure 3(c)), and will merge with the 

290 MCl-driven shock into a stronger compound shock, which is consistent with the previous 

291 results of double MC interaction [Odstrcil et ai, 2003; Lugaz et ai, 2005]. Therefore the 

292 Multi-MC has nearly been completing its final evolutionary stage at 53.3 hours, after 

293 which the Multi-MC will move forwards as a relatively stable structure. 

294 Time sequence of hypothetical measurement at LI for Case Ci is shown in Figure 4. 

295 The MC2-driven shock just emerges from MCI body after penetrating it, so no extremum 

296 of speed profile Vr is found inside the multi-cloud. Double dips of Dst index are —93 nT 

297 and —95 nT, increased by 47 nT and decreased by 20 nT respectively in contrast with 
29E those in Case Bi in Figure 2(e). The mitigation of geoeffectiveness for the first Dst dip 

299 is owing to the position of MC2-driven shock front far away from the rear part of MCI 

300 with southward magnetic component, the aggravation for the second Dst dip is ascribed 

301 to the MC2 body compression mentioned above. A peak of VB^ up to 14 mV/m can be 

302 seen near the MCI front boundary, where the largest compression occurs. However, it 

303 is positive and makes no contribution to geoeffectiveness. Additionally, the durations of 
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304 MCI and MC2 are shortened by 4.7 and 3 hours, respectively, as compared with those in 

305 Case Bi. 

306 Figure 5 shows the time-dependent parameters of Multi-MC Case Ci (thick curves), 

307 where the dotted, dashed, and dotted vertical lines from left to right denote the occasion 

308 of MC2-driven shock encountering MCI body tail, MC2 body hitting MCI body tail, and 

309 MC2-driven shock reaching MCI body head, respectively. Two corresponding isolated 

310 MC cases are superimposed as thin curves for comparison. The acceleration of MCI is 

311 large and early, while the deceleration of MC2 is small and late, as seen from Figure 5(a). 

312 The radial compression of MC2 body brings not only the shortening of its radial span 

313 Sr but also the stretching of its angular span 5*^. The behavior of MCI is a bit more 

314 complex. In our previous studies of MC-shock interaction [Xiong et ai, 2006a, b], MC 

315 compressed morphology will be restored after the shock passage. However, in the presence 

316 of the following MC2 body's pushing effect for Multi-MC case, MC2 body will take over 

317 the role of suppressing MCl's inherent magnetic elasticity when MC2-driven shock moves 

318 farther. Thus the firm gripping of MCI body at all time leads to significant shrinking of 

319 its volume. It is why cross section area of MCI body is smaller than that of MC2 (Figure 

320 5(d)). Particularly, as seen from the local minimum value of Sr at 30 hours in Figure 

321 5(b), the compression of MCl's Sr reaches to its extreme when the MC2-driven shock 

322 nearly arrives at MCI head boundary. Meanwhile the temporarily enhanced SO of MCI 

323 during 24 ~ 38 hours is steadily reduced afterwards (Figure 5(c)). 

3.3. MC2 Helicity Role 

324 There are various combination modes to form a double-MC structure on basis of each 

325 sub-cloud helicity signature [Wang et ai, 2002], one of which possessing the strongest 
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326 geoeffectiveness is positive helicity for preceding MCI (i/mci = 1) and negative helicity 

327 for following MC2 {Hmc2 = —1) [Wang et ai, 2004]. According to this scenario [Wang 

328 et ai, 2004], simulation cases B2 and C2 are run simply by reversing MC2 helicity in Cases 

329 Bi and Ci, respectively. The in-situ record of passage of multi-cloud event at LI is shown 

330 in Figure 6, with Columns (A) and (B) corresponding to Cases B2 and C2, respectively. 

331 In contrast to Figures 2 and 4, the elevation angle 6 of magnetic field vector within 

332 the double-fiux-rope structure in Figure 6 is changed from the north-south-north-south 

333 orientation to north-south-south-north one. Though two Dst dips exist in Groups EID2 

334 and CID2, close scrutiny reveals that (1) the recovery phase of the first trivial Dst dip 

335 is extremely short (3.3 and 0.9 hours in Cases B2 and C2, respectively); (2) the second 

336 Dst dip is low enough to describe the whole geoeffectiveness by its local minimum, with 

337 —166 nT at 90 hours in Case B2 and —144 nT at 78 hours in Case C2. Hence from the 

338 perspective of continuous interval with southward magnetic field Bg, Dst curve in Groups 

339 EID2 and CID2 can be considered as a one-dip structure by ignoring the first trivial 

340 dip. The closer the distance between two sources of IP geoeffective trigger, the easier is 

341 the superposition of individual geoeffectiveness, the greater is the resulting geomagnetic 

342 storm. This is confirmed by contrast of Figures 6(c) and (f) with Figures 2(e) and 4(e). 

4. Geoeffectiveness Studies 

343 Near-HCS latitudinal dependence of the Dst index is plotted in Figure 7, where Columns 

344 (A) Dstpi and (B) Dstp2 represent the first and the second Dst dips in Cases Bi and 

345 Ci, meanwhile Column (C) DstN depicts the single Dst dip in Cases B2 and C2. The 

346 dashed and dash-dotted lines represent for Cases Bi and Ci, respectively, in Columns 

347 (A) and (B). And they represent for Cases B2 and C2, respectively, in Column (C). As 
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348 the MC2-driven shock continues to propagate through the MCI medium, Dstpi increases 

349 within Lat. > 1.3° and decreases within Lat. < 1.3°, found in Figure 7(A). Meanwhile 

350 the distribution of Dst^ in Figure 7(C) is quite similar. The trend of decreased Dst 

351 near HCS is opposite to that in the case of MC-shock interaction [c.f. Figure 8 in Xiong 

352 et al, 2006a]. The above divergence is clarified by the absence of following MC body 

353 pushing in MC-shock interaction [Xiong et al., 2006a]. First, the latitudinal extent of 

354 MC2 body is much narrower than that of MC2-driven shock. Second, the coalescent 

355 boundary between MCI and MC2 body is further narrower, which covers latitude range 

356 between 4.5°S and 4.5°N (Figure 3(c)). Thus the MC2 body pushing effect is strongest at 

357 the equator, within confined latitudinal extent between 4.5°S and 4.5°N. The near-HCS 

358 geoeffectiveness of Dstpi from Case Bi to Case Ci, DsIn from Case B2 to Case C2 is 

359 subsequently aggravated. As a result, nonuniform latitudinal distribution of Dstpi and 

360 Dst]sr is intensified. Besides, Dstp2 is nearly unaffected in Case Bi. However Dstp2 in 

361 Case Ci is obviously decreased, as a result from the compression of MC2 body interpreted 

362 in Section 3.2. Hence the geoeffectiveness of Multi-MC is indeed largely enhanced due to 

363 interaction between sub-clouds, as compared with that in an isolated MC event. 

364 In order to quantify the evolution process of Multi-Cloud, do = r^ci — '^'mcii the distance 

365 between the cores of MC2 and MCI on the occasion of MCI head just reaching LI, is 

366 chosen as an indicative parameter. Tmci and rmc2 are the core positions of MCI and MC2 

367 in radial direction respectively. The reliance of several multi-cloud parameters on do is 

368 further explored in Figure 8 by the integrated study of Groups EIDi and EID2. The 

369 absolute value of is labeled as |do|- As Dt, the emergence interval of MCI and MC2, 

370 decreases, |(io| firstly reduces from 107 to 53 Rg at a constant slope, then asymptotically 
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371 approaches to 42 Rs shortly after MC2-driven shock emerges from MCI body head (Figure 

372 8(a)). The penetration depth of MC2-driven shock in MCI medium, dj^gt, defined by 

373 the radial distance between MC2-driven shock front and MCI inner boundary along the 

374 equator, is shown in Figure 8(b), which can be divided into four stages according to 

375 the different behaviors: (1) a rapid increase during do < — 66.3/?^, (2) an extremely 

376 slow increase during — 66.3i?s < do < — 52.5i?s, (3) a fast re-increase during —52.5Rs < 

377 do < — 4:6.7 Rs, (4) a very small oscillation around the final limit value of 40 Rg during 

378 do > — 46.7i?s- The rapid increasing of dost in stages (1) and (3) is straightforward due 

379 to continuous forward movement of shock front in MCI medium. As for stage (2) during 

380 which the shock front hits MCI core, though the shock front location relative to the 

381 MCI body is deeper and deeper at that time, the abrupt change of MCI rear boundary 

382 morphology from a V-shape to a straight line, mentioned in Section 3.2, greatly reduces the 

383 radial extent of MCI rear half, and hence significantly inhibits the increase of the absolute 

384 value of penetration depth d^st- When the shock front crosses the MCI front boundary 

385 {do > — 46.7 Rs), the magnetic tension of the highly compressed MCI body is drastically 

386 accumulated. As a result, the nearer the distance between two sub-MCs is (the shorter 

387 the I do I is), the larger is the resistance of MCI elasticity against compression. The final 

388 equilibrium is naturally manifested in the behavior of stage (4). The early and sensitive 

389 response of Max.(i?mci) at do = — lOOi?^ is conspicuous along Lat. = 4.5° in Figure 8(c), 

390 because the initial interaction between MCI body and MC2-driven shock happens around 

391 Lat. = 4.5°. The swift enhancement of Max.(i?mci) during — 68.9/?^ < do < —46.7Rs 

392 along Lat. = 0° is owing to the compression concurrently exerted by the MC2-driven 
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393 shock and MC2 body. Both Max.(i?mci) and Max. (5^02) reach a relatively stable state 

394 when Idol = 42i?<.. 

395 The variance of geoeffectiveness as a function of rfo is elucidated in Figure 9. The 

396 analyses on Dstpi, Dstp2, and DstN are addressed one by one. First, when do < —GORg, 

397 the behavior of all parameters in Figures 9(a)-(d) for Dstpi is pretty coincident with 

398 that of our previous study for MC-shock interaction [c.f. Figures 9(b)-(e) in Xiong et al, 

399 2006a]. The dynamics of MC1-MC2 merging at that time is dominated by the interaction 

400 between MC2-driven shock and MCI body. Thus MC2-driven shock plays the similar 

401 role of the incidental shock as addressed before [Xiong et ai, 2006a], which clarifies the 

402 above-mentioned coincidence. As \do\ is reduced from 60Rs to 52.5Rs, MC2 body directly 

403 collides with MCI body. It leads to the decrease of Dst, Min.{VBz), and Min.(i?s) due 

404 to compression. Particularly, the decrease of Min. (VBz) and Mm.{Bs) along Lat. = 

405 4.5°S is very drastic, because the change of MCI field line morphology from a V-shape 

406 to a straight line mentioned in Section 3.2 leads to the southward rotation of magnetic 

407 field within MCI rear half along Lat. = 4.5°S. This additional rotation effect further 

408 strengthens Bg along Lat. = 4.5°S. When |(io| continues to decrease to be less than 52.5/2^, 

409 significant difference of geoeffectiveness between Lat. = 0° and 4.5°S occurs. Along Lat. 

410 = 4.5°S, the rapid recovery of Min.(Ss) from —24.5 to —13.5 nT, and Mm.{VBz) from 

411 —15 to —8 mV/m, leads to the subdued Dstpi from its minimum —165 to —100 nT. 

412 Contrarily, the geoeffectiveness along Lat. = 0° remain unchanged (Figures 9(a)-(d)). 

413 Namely, the aggravated geoeffectiveness along the equator is the same with Dstpi = 

414 —180 nT, provided that \do\ is smaller than a certain threshold of 52.5 Rg. This highly 

415 nonuniform latitudinal distribution of Dstpi is owing to the limited latitudinal range 
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416 (4.5°S ~ 4.5°N) of pushing effect of MC2 body. When the shock ultimately penetrates 

417 MCI body, the persistent pushing of following MC2 body within 4.5°S ~ 4.5°N can prevent 

418 the previously compressed magnetic field lines of MCI body from being relaxed. So Dstpi 

419 along Lat. = 0° is nearly constant for \do\ < 52.5i?s. As for Dstpi along Lat. = 4.5°S, it 

420 increases as a result of relaxation of magnetic tension without MC2 body pushing. Second, 

421 the variance of Dstp2 (geoeffectiveness of sub-MC2) only happens between do = — 68i?s 

422 and — 46.7i?s, during which the MC2 body compression due to the blocking of MCI 

423 body takes effect. Before the involving of MC2 body into interaction {do < —68Rs), 

424 or after the completion of Multi-MC's drastic evolution stage {do > — 46.7/?^), Dstp2 

425 is unchanged. By comparison Dstp2 with Dstpi, one can see that the MCI undergoes 

426 the greater compression than the MC2. Third, the behavior of DsIn (Figures 9(i)-(l)) 

427 is quite similar to that of Dstpi (Figures 9(a)-(d)) due to similar reasons mentioned 

428 above. The minimum Dst in Figures 9(a), (e), and (i) is —180, —130, and —235 nT, 

429 respectively. The greatest geoeffectiveness of Dstj\f directly results from the longest At 

430 (Figure 9(k)). Therefore, the geoeffective parameters of every sub-MC are dramatically 

431 changed in contrast with those of the corresponding isolated MC during the merging 

432 process. For the IP compound structure formed by multiple ICMEs, the geoeffectiveness 

433 is jointly determined by two factors: the parameters of the individual ICMEs themselves, 

434 and the interaction process between these ICMEs. This is substantiated by the observation 

435 data analyses [Wang et al, 2002, 2003a; Xue et al, 2005; Farrugia et al, 2006; Zhang 

436 et ai, 2007] and our quantitative investigation of numerical simulation of this study. 

437 The Multi-Cloud geoeffectiveness depends on not only the MC1-MC2 eruption interval, 

438 but also collision intensity. Obviously, an MCI overtaken by an MC2 with various initial 
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439 speeds may result in different geoeffectiveness. From the Figure 9 concerning Groups EIDi 

440 and EID2, two basic results are obtained: (1) The maximum geoeffectiveness occurs at Lat. 

441 = 0° for the same propagation direction of MCI and MC2 along the equator; (2) The final 

442 Dst at Lat. = 0° is nearly constant, provided the accompanying |(io| is sufficiently small 

443 (|(io| < 46.7-Rs), or the initial MC1-MC2 eruption interval is sufficiently short {Dt < 20 

444 hours). With tmc2 designated to be 12.2 hours, the reliance of geoeffectiveness along the 

445 equator on collision degree is further explored in Figure 10 by parametric study of variable 
«6 Vmc2- The larger the value of Vmc2 is, the greater is the collision degree that the Multi-MC 

447 may suffer from. Dstpi only decreases a bit from —180 to —210 nT within such a wide 

448 spectrum of Vmc2 from 450 to 1200 km/s. The geoeffectiveness enhancement of Multi- 

449 MC is ascribed to compression between the sub-MCs. When the MCI compression has 

450 already approached to saturation, the effect to increase MCI geoeffectiveness by having 

451 MCI impinged by a highly fast MC2 is extremely limited. It is more and more difficult to 

452 quench the dramatically accumulated magnetic elasticity of MCI body, as MCI undergoes 

453 the greater and greater compression. The impact of the high-speed MC2 body is largely 

454 offset by the buffer action of magnetic tension of the MCI body. As for Dstp2, the increase 

455 of Vmc2 has a direct influence. However, Dstp2 deducted by the Dst of the corresponding 

456 individual MC2 event is roughly constant, which can be seen from Figure 10(b). Namely 

457 Dstp2 decreases from —125 to —190 nT, as Vmc2 increases from 450 to 1200 km/s, chiefly 

458 ascribed to the increase of geoeffectiveness of the corresponding individual MC2 event 

459 itself, but not MC1-MC2 interaction. Excluding the geoeffectiveness increase of individual 

460 MC2 event, Dst^ still decreases for Vmc2 > 1000 km/s in Figure 10(c), because interaction 

461 takes obvious effect herein. The geoeffectiveness variance can be elucidated from the 
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462 perspective of dynamic response of sub-MCs. The double-MC interacting region is within 

463 MCI rear part and MC2 front part, where the direct compression occurs. So the factor of 

464 MC1-MC2 interaction for geomagnetic storm enhancement is strongest for Dstj^i, weakest 

465 for Dstp2- In conclusion, two points can be drawn from Figures 9 and 10: (1) The 

466 significant geoeffectiveness variance accompanies the different evolution stages; (2) Once a 

467 Multi-MC completes its evolution process before its arrival at 1 AU, the collision intensity 

468 between sub-MCs merely modulates the final geoeffectiveness a bit. The innate magnetic 

469 elasticity can buffer the reciprocal collision between sub-MCs against each other. When 

470 every sub-MC becomes stiffer and stiffer, the compression reaches its asymptotic degree, 

471 and the geoeffectiveness enhancement becomes less and less obvious. Therefore, with 

472 respect to Multi-MC geoeffectiveness, the evolution stage is a dominant factor, whereas 

473 the collision intensity is a subordinate one. 

474 Additionally, the dependence of geoeffectiveness of an individual MC on the eruption 

475 speed Vmc is also revealed from the isolated MC2 events from Figure 10. If region in 

476 MC medium is located in its anterior half (Group IM2), Dst steadily decreases as Vmc 

477 increases, as seen by the thin solid line in Figure 10(C); Contrarily, if Bg region is to be 

478 in the rear half of MC (Group IMi), Dst only decreases on the condition of Vmc > 800 

479 km/s, as seen by the solid line in Figure 10(B). The increase of Vmc leads to a more violent 

480 interaction of individual MC body with the ambient solar wind ahead. As a result, MC 

481 core, initially located at the geometry center of MC boundary, will be gradually shifted to 

482 MC anterior boundary. MC anterior half is preferential compressed, because MC-ambient 

483 flow interaction originates from MC front boundary. The compression exists in MC rear 

484 half, only when the whole cross section area of MC body is significantly contracted on the 
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485 condition of very fast speed Vmc- This is why Dst for Group IMi remains a constant of 

486 -100 nT within = 450 ~ 800 km/s. 

5. Compressibility Analyses 

487 The idea that the compression is an efficient mechanism to enhance the geoeffective- 

488 ness of the pre-existing Bg event has been proved in data analyses [Wang et ai, 2003c]. 

489 Compression effect is virtually responsible for the geoeffective property of Multi-MC. So 

490 it is very meaningful to analyze the maximum compression degree for a Multi-MC. 

491 The Multi-MC characteristics can be inferred from several parameters of near-Earth 

492 measurements, depicted by Figure 11. The interchange of momentum between the pre- 

493 ceding slow cloud MCI and following fast cloud MC2 leads to MCI acceleration and MC2 

494 deceleration, which influences, more or less, the Sun-Earth transient time, TT^ci and 

495 TTmc2 for MCI and MC2, respectively. The shortening of TTmci begins at Dt = 21 hours 

496 as seen from Figure 11(a), meanwhile the lengthening of TTmc2 begins at Dt = 28 hours, 

497 seen from Figure 11(b). The MCI acceleration is very obvious, as the larger Vmc2 is? the 

498 smaller is TT^ci (Figure 11(g)). Contrarily, the MC2 slowdown is independent of Vmc2, 

499 as TTmc2 in coupled events deviates from that in the corresponding isolated events by a 

500 constant (Figure 11(h)). The effect oiTTmci decrease is much greater than that of TTmc2 

501 increase for Multi-MC cases. Since the transporting time of an ICME may be modified if 

502 it interacts with others during its IP propagation, some empirical formulas of transporting 

503 time on basis of observations of one single ejecta event [Gopalswamy et ai, 2000, 2001a] can 

504 not be directly applied to the ICME-ICME interaction cases [Farrugia and Berdichevsky, 

505 2004; Wang et ai, 2005b; Xiong et ai, 2005]. Coupling between ICMEs occupies a large 

506 fraction for the causes of great geomagnetic storms [Xue et ai, 2005; Zhang et ai, 2007], 
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507 the Multi-MC should be paid special attention for space weather predicting. Thus the 

508 numerical simulation based on physics models is very useful to forecast the arrival time 

509 of the interacting ICMEs. The duration of sub-MC passage at LI, AT^c, is a distinct 

510 reflection of compression effect. AT^ci exists a lower limit, as shown in Figure 11(c), so 

511 does ATmc2 in Figure 11(d). When the Multi-MC experiences the sufficient evolution for 

512 Dt < 24 hours, the reduction of ATmci and ATmc2 is 14 and 4.5 hours respectively, in 

513 contrast with the corresponding isolated sub-MC cases. As Vmc2 increases, both ATmci 

514 and ATmc2 monotonically decrease. However, the solid and dashed lines, representing the 

515 Multi-MC and corresponding isolated MC events in Figure ll(j), intersect at 'i'm,c2 — 

1040 

516 km/s. ATmc2 in MC1-MC2 interaction is determined by two factors: (1) the compression 

517 of MC2 radial extent resulting from collision; (2) the slowdown of MC2 body as a result of 

518 momentum transfer from MC2 to MCI body. The first factor, tending to shorten ATmc2, 

519 dominates the cases for Vmc2 < 1040 km/s; the second factor, trending to lengthen ATmc2, 

520 dominates the cases for Vmc2 > 1040 km/s. Besides, the near-Earth radial span of MCI 

521 body Srmci in Figures 11(e) and (k) has the similar variance trend as ATmci in Figures 

522 1 1 (c) and (i). It again proves that the compression has saturation effect for MCI body. 

523 The Srjnci of 67 Rg in an individual case can be compressed to 40 Rg at most by Dt 

524 reduction (Figure 11(e)). Sv^ci decreases very slowly from 43 Rg at Vmc2 ~ 

450 km/s 

525 to 25 Rg at Vmc2 = 1200 km/s (Figure ll(k)). Moreover, the overall compression degree 

526 for a Multi-MC is well described by do, the distance between the core of the following 

527 MC2 and preceding MCI on the occasion of MCI head just reaching LI. One can see 

528 that do variance is associated with Multi-MC evolution stages (Figure 11(f)). The swiftly 

529 reducing trend of \do\ at the beginning is suddenly stopped at Dt = 25 hours. \do\ reaches 
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530 its lower limit of 42 Rs at Dt =17 hours, and maintains a horizontal slope afterwards. 

531 When the inherent magnetic tension rivals the external compression for force balance, 

532 each sub-MC behaves like a rigid body with a little elasticity. \do\ is only reduced from 

533 43 to 30 Rs over such a wide Vmc2 range from 450 to 1200 km/s (Figure 11(1)). 

534 The compression due to interaction is primarily responsible for geoeffectiveness en- 

535 hancement, once two MCs form a Multi-MC. Assuming nonexistence of magnetic field 

536 in the IP medium and all ejecta, the preceding ejecta may be exorbitantly compressed 

537 to an unbelievably small scale by the following ejecta [Gonzalez- Esparza et al, 2004; 

538 Gonzalez-Esparza, 2005]. Obviously the compressibility on basis of hydrodynamic nature 

539 [Gonzalez-Esparza et al, 2004; Gonzalez-Esparza, 2005] is overestimated due to ignoring 

540 of magnetic elasticity. The larger the compression is, the stiffer is every sub-MC body. 

541 Hence a cutoff compression degree exists because of magnetic tension. Besides, if the 

542 helicity of MCI is consistent with that of MC2, a electric current sheet occurs between 

543 the adjoining boundary of MCI and MC2 due to magnetic field direction reversion. The 

544 electric current intensity synchronously increases with the Multi-MC compression. If mag- 

545 netic reconnection happens there, the MC1-MC2 collision effect would be weakened. As 

546 a result, the outermost part of magnetic field lines of each sub-MC would be reconnected 

547 together [Wang et al., 2005b]. Particularly in the condition of large speed difference be- 

548 tween MCI and MC2, both MCs may be merged into one new magnetic flux rope by 

549 the driven magnetic reconnection [Odstrcil et al, 2003; Schmidt and Gargill, 2004; Wang 

550 et al, 2005a]. The magnetic reconnection reduces Multi-MC's cutoff compression degree. 

551 If magnetic reconnection is introduced into Groups EIDi and CIDi of Table 2, the Multi- 

552 MC geoeffectiveness would become weakening due to the subdued compression and south 
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553 magnetic component annihilation. However, magnetic diffusion in the IP space should be 

554 very small, magnetic reconnection may sightly modulate, but not significantly distort the 

555 dynamics and geoeffectiveness of Multi-MC in the framework of ideal MHD process. So 

556 the CME-CME cannibalization, firstly observed in the inner corona by the SOHO/Lasco 

557 [Gopalswamy et ai, 2001b], later proved to be caused by magnetic reconnection [Wang 

558 et ai, 2005a], may not occur in the IP space [Wang et ai, 2005b]. 

6. Conclusions and Summary 

559 In order to better understand the nature of IP Multi-MC structure, the interaction 

560 between two IP MCs (MCI and MC2), and the ensuing geoeffectiveness are explored 

561 under a very simplified and specialized circumstance by a 2.5-dimensional ideal MHD 

562 numerical model. This work is a continuation to our recent studies of MC-shock interaction 

563 [Xiong et ai, 2006a, b] by replacing a following incidental strong shock with a following 

564 fast MC. Via analyses of a comprehensive integration of many simulation cases under 

565 various conditions, it is found that the magnetic elasticity, magnetic helicity of each MC, 

566 and compression between each other are the overriding physical factors in the formation, 

567 propagation, evolution, and resulting geoeffectiveness of IP Multi-MC. 

568 First, the dynamical response of MCs colliding is studied. The coupling of two MCs 

569 could be considered as the comprehensive interaction between two systems, each com- 

570 prising of an MC body and its driven shock. Because the following MC2 is faster than 

571 the preceding MCI, the MC2-driven shock and MC2 body successively impact the rear 

572 boundary of MCI body. As a result, the morphology of magnetic field lines at MCl's 

573 rear part is consequently changed from its initial rough semi-circle to a V-shape, and 

574 then to a flat line. As swept by the marching MC2-driven shock front, the local mag- 
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575 netic field lines in MCI medium just downstream of MC2-driven shock front would be 

576 compressed and rotated. The pushing of MC2 body prevents the previously compressed 

577 magnetic field in MCI medium from being restored, after the passage of MC2-driven shock 

578 front. MCI body undergoes the most violent compression from the ambient solar wind 

579 ahead, continuous penetration of MC2-driven shock through MCI body, and persistent 

580 pushing of MC2 body at MCI tail boundary, which leads to a significant shrinking of 

581 MCl's cross section. Contrarily, the blocking of MCI body also results in the change of 

582 MC2 boundary from a radial-extent-elongated ellipse to an angular-extent-elongated one. 

583 The Momentum is continuously transferred from sub-MC2 to sub-MCl, until the radial 

584 profile of Multi-MC speed is monotonically decreasing with the maximum value at MCl- 

585 driven sheath. When MC2-driven shock is merged with MCl-driven shock into a stronger 

586 compound shock, Multi-MC completes its ultimate evolutionary stage, and hence moves 

587 forward relatively stable entity. 

588 Second, the geoeffectiveness of MCs coupling is explored. The interaction of MCI 

589 and MC2 in the IP space results in the superposing of their geoeffectiveness. The two- 

590 MC event is associated with a two-step geomagnetic storm, as indicated by two Dst dips. 

591 Particularly, if Bg region in a Multi-MC is located at MCI rear half and MC2 anterior half, 

592 the Multi-MC excites the greatest geomagnetic storm among all combinations of each sub- 

593 MC helicity, and two Dst dips can be nearly reduced to a single Dst dip due to ignoring of 

594 the very short recovery phase of the first Dst dip. The geoeffectiveness of each individual 

595 MC is largely enhanced as a result of MC1-MC2 interaction. Moreover, because latitudinal 

596 extent of MC body is much narrower than that of its driven shock, the effect of MC2 body 

597 pushing upon MCI body is limited within a very narrow latitudinal band centered at the 
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598 heliospheric equator. Outside this latitudinal band, geoeffectiveness is initially enhanced 

599 and then recovered, as the emergence interval of two MCs becomes shorter and shorter; 

600 meanwhile the geoeffectiveness is firstly aggravated and then maintains constant inside 

601 this band. Obviously, the nonuniform latitudinal distribution of geoeffectiveness is further 

602 intensified by MC2 body pushing. Moreover, With respect to Multi-MC geoeffectiveness, 

603 the evolution stage is a dominant factor, whereas the collision intensity is a subordinate 

604 one. 

605 Third, Multi-MC's compressibility associated with magnetic elasticity is analyzed. Both 

606 compression degree and evolutionary stage of a Multi-MC could be quantitatively de- 

607 scribed by \do\, the absolute distance between MCI and MC2 core on the occasion of 

608 MCI head just reaching LI. The shorter the \do\ is, the greater is Multi-MC's compress- 

609 ibility. Magnetic field lines of MCI body initially appears to be too frail to resist the 

610 collision in the face of the overtaking high-speed MC2, so Idol is steadily reduced. As the 

611 evolution of Multi-MC proceeds, the MCI body suffers from larger and larger compres- 

612 sion, and its original vulnerable magnetic elasticity becomes stiffer and stiffer. When the 

613 accumulated inherent magnetic elasticity in the highly shrunk MCI body can counteract 

614 the external compression, the previous continuously reducing |(io| drastically approximates 

615 to an asymptotic limit. Magnetic elasticity not only buffers the collision between MCs, 

616 but also leads to a cutoff compression degree of Multi-MC. Moreover, the collision of MC2 

617 with a very wide speed spectrum upon MCI has a little influence to enhance the cutoff 

618 compressibility. However, if magnetic reconnection occurs within the interacting region 

619 of Multi-MC, the cutoff compressibility would be expected to decrease a bit. 
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620 Overall, the Multi-MC is of great concern for space weather community. The geoef- 

621 fectiveness enhancement of coupling of multiple MCs is virtually ascribed to compression 

622 in the Multi-MC. The maximum compressibility of Multi-MC is mainly decided by its 

623 inherent magnetic elasticity. 
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Figure Captions 

874 Figure 1 The evolution of MC2 overtaking MCI for Case Bi, with (a)-(c) magnetic 

875 field magnitude B, (d)-(f) radial flow speed Vr, and (g)-(i) radial characteristic speed of 

876 fast mode c/. Attached below each image are two additional radial profiles along Lat.= 0° 

877 and 4.5°S. Note that radial profile of B is plotted by subtracting the initial ambient value 

878 B\t=o. The white solid line in each image denotes the MC boundary. Solid and dashed 

879 lines at each profile denote MC core and boundary. Only part of domain is adaptively 

880 plotted to highlight Multi-MC. 

881 Figure 2 In-situ hypothetical observation along Lat. = 4.5°S for Case Bi. Stacked 

882 from top to bottom are (a) magnetic field magnitude B, (b) elevation of magnetic field 

883 G, (c) radial flow speed Vr, (d) derived dawn-dusk electric field VB^, and (e) Dst index. 

884 Solid and dashed delimiting hues denote MC center and boundary. 
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S85 Figure 3 The evolution of MC2 overtaking MCI for Case Ci with radial flow speed Vr- 

886 Figure 4 In-situ hypothetical observation along Lat. = 4.5°S for Case Ci. 

887 Figure 5 The time dependence of MC parameters: (a) radial distance of MC core r^, 

888 (b) MC radial span Sr, (c) MC angular span S9, and (d) MC cross section area A. The 

889 thick dashed and solid lines denoted the preceding MCI and following MC2 in Multi- 

890 MC Case Ci, superimposed with thin lines for corresponding individual MC cases for 

891 comparison. Three vertical delimiting lines (dotted, dashed and dotted) from left to right 

892 correspond to the occasion of MC2-driven shock encountering MCI body tail, MC2 body 

893 hitting MCI body tail, and MC2-driven shock reaching MCI body head, respectively. 

894 Figure 6 In-situ hypothetical observation along Lat. = 4.5°S for (A) Case B2 and (B) 

895 Case C2. Cases B2 and C2 differ from their respective companion Cases Bi and Ci by the 

896 opposite MC2 magnetic helicities. 

897 Figure 7 The comparison of latitudinal distribution of Dst index among the Multi-MC 

898 Cases Bi, Ci, B2, and C2. Double Dst dips in Cases Bi and Ci with positive magnetic 

899 helicities in MC2 are shown in (A) Dstpi and (B) Dstp2, as well as a single Dst dip 

900 in Cases B2 and C2 with negative helicity in MC2 shown in (C) Dstj^. Dashed and 

901 dash-dotted lines in (A) and (B) correspond to Cases Bi and Ci respectively; dashed and 

902 dash-dotted lines in (C) correspond to Cases B2 and C2 respectively. The isolated events 

903 corresponding to MCI and MC2 for Case Bi are denoted as solid lines in (A) and (B), 

904 those for Case B2 as solid thick and thin lines in (C). 

905 Figure 8 The do-dependent parameter variances at LI in Group EIDi: (a) Dt, time 

906 interval of MCs {Dt = tm.c2 — tmci, tmci = hour); (b) dost, penetration depth of MC2- 

907 driven shock in MCI medium; (c) Max.{Bmci), the maximum of magnetic field strength 
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908 in MCI; and (d) Max.(i?mc2), the maximum of magnetic field strength in MC2. Here 

909 do refers to the distance between MC2 core rmc2 and MCI core r^ci on the occasion of 

910 MCI head just reaching LI, namely do = rmc2 ~ fmci- The vertical delimiting dotted and 

911 dashed lines denote the occasions of MC2-driven shock just hitting MCI core and head 

912 at LI. In (c) and (d), the thick solid and dashed lines denote observations along Lat. = 

913 0° and 4.5°S, while the thin horizontal ones represent the values of corresponding isolated 

914 MC events. 

915 Figure 9 The parameter variances of Multi-MC geoeffectiveness as a function of d^. (a, 

916 e, i) Dst index; (b, f, j) Min.(yi?^), the minimum of dawn-dusk electric field VBz'., (c, g, k) 

917 At, the interval between the commencement of VBz < —0.5 mV/m and the corresponding 

918 minimum Dst; and (d, h, 1) Min.(i?s), the minimum of southward magnetic component. 

919 Solid and dashed lines correspond to observations along Lat.= 0° and 4.5°S respectively. 

920 The double Dst dips in Group EIDi are shown by Columns (A) Dstpi and (B) Dstp2, 

921 and the single Dst dip in Group EID2 by Column (C) Dst^- The horizontal solid and 

922 dashed lines denote observations of the isolated events, corresponding to Group EIDi, at 

923 Lat. = 0° and 4.5°S respectively, with MCI in Column (A) and MC2 in Column (B). 

924 Figure 10 The reliance of Dst in Multi-MC cases on initial speed of following MC2. 

925 Double Dst dips in Group CIDi are shown as dashed lines by (A) Dstpi and (B) Dstp2, 

926 while a single Dst dip in Group CID2 by (C) Dst^- Decoupled MCI and MC2 events in 

927 Group CIDi are plotted as solid lines in (A) and (B), those in Group CID2 as thick and 

928 thin solid lines in (C). The curves of single MC2 event in (B) and (C) are non- horizontal 

929 due to Vmc2 variance. 
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930 Figure 11 The dependence of Multi-MC characteristic parameters at LI on MC2-MC1 

931 eruption interval —Dt {—Dt = — 1 ■ Dt = tmci — tmc2) and MC2 speed Wmc2 is shown by 

932 Columns (A) Group EIDi and (B) Group CIDi. TT^ci, the Sun-Earth transient time of 

933 MCI, is shown in (a) and (g); AT^ci, MCI event duration at LI, in (c) and (i); S'r.mci, 

934 MCI radial span, in (e) and (k); and d^, the distance between MC2 and MCI core, in 

935 (f) and (1). (a), (c), (g), (i) and (b), (d), (h), (j) are the counterparts for MCI and MC2 

936 respectively, (e), (f), (k), and (1) refer to the occasion when MCI head just reaches LI. 

937 Dashed lines in all panels except (f) and (1) represent the corresponding isolated MC 

938 events for comparison. The vertical dotted and dashed lines in Column (A) denote the 

939 cases of MC2-driven shock just hitting MCI core and head respectively. 
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Table 1. Assortment of simulation cases of individual MC 



Group 


Case 


(10^ km/s) 


Comment 


IMi 


bi, ci, di, 


4, 6, 5, 


Individual MC 




ei, fi, gi, 


7, 8, 9, 


{Hmc = 1) 




hi, ii, ji 


10, 11, 12 




IM2 


b2, C2, d2, 


4, 6, 5, 


Individual MC 




e2, f2, g2. 


7, 8, 9, 


(Hmc = — 1) 




h2, i2, h 


10, 11, 12 
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Table 2. Assortment of simulation cases of Multi-MC. Note that Vmci = 400 km/s, 



tmci = hour for all 48 cases. 



Group 


Case 


Vmc2 

(10^ km/s) 


tmc2 

(hour) 


Comment 


ElDi 


Bi, Ci, Di, El, 
Fi, Gi, Hi, li, 
Ji, Ki, Li, Ml, 
Ni,Oi, Pi,Qi 


6 


30.1, 12.2, 44.1, 42.1, 
40.1, 37.1, 35.1, 33.1, 
31.5, 28.1, 25.1, 22.1, 
20.1, 17.1, 15.1, 10.2, 


Eruption Interval 
Dependence 

{Hmcl = 1, Hmc2 = 1) 


EID2 


B2, C2, D2, E2, 
F2, G2, H2, I2, 
J2, K2, L2, M2, 
N2, O2, P2, Q2 


6 


30.1, 12.2, 44.1, 42.1, 
40.1, 37.1, 35.1, 33.1, 
31.5, 28.1, 25.1, 22.1, 
20.1, 17.1, 15.1, 10.2, 


Eruption Interval 
Dependence 

{Hmcl = 1; Hmc2 = ~1) 


ClDi 


Ri, Si, Ci, Ti, 
Ui, Vi, Wi, Xi, 
Yi 


4.5, 5, 6, 7, 
8, 9, 10, 11, 
12 


12.2 


Collision Intensity 
Dependence 

{Hmcl = I7 Hmc2 = 1) 


CID2 


R2, S2, C2, T2, 
U2, V2, W2, X2, 
Y2 


4.5, 5, 6, 7, 
8, 9, 10, 11, 
12 


12.2 


Collision Intensity 
Dependence 

{Hmcl = 1; Hmc2 = ^1) 
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Figure 2. 
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Figure 3. 
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In-Situ Observation at LI along Lat. = 4.5°S 




Time (Hour) 



Figure 4. 
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Figure 6. 
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Figure 7. 
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Figure 9. 
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Figure 10. 



DRAFT 



April 5, 2009, 4:06ain 



DRAFT 




Figure 11. 



DRAFT 



April 5, 2009, 4:06ain 



DRAFT 



